IC 4406: a radio- infrared view 
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new observations of this source performed with IRAC (Spitzer Space Telescope) 
and the Australia Telescope Compact Array. The radio maps show that the flux 
from the ionized gas is concentrated in the bright central region and originates 
l/^ | in a clumpy structure previously observed in Ha, while in the infrared images 

filaments and clumps can be seen in the extended nebular envelope, the central 
region showing toroidal emission. Modeling of the infrared emission leads to the 
conclusion that several dust components are present in the nebula. 
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ABSTRACT 

IC 4406 is a large (about 100" x 30") southern bipolar planetary nebula, 
composed of two elongated lobes extending from a bright central region, where 
there is evidence for the presence of a large torus of gas and dust. We show 



Subject headings: Radio continuum: stars, Infrared: stars, Planetary nebulae: 
individual (IC 4406) 



1. Introduction 

IC 4406 is a well-studied southern planetary nebula. It has been image d with severa l 



telescopes at different wavelength ranges. Near-IR images show two H 2 lobes (jStoreylll984l ) 



orthogonal to the nebula's major axis and ~25" away from each other. Thes e peaks are 



appro ximately coincident with the two blobs observed in Ha+[N II] and [O III] ( iSahai et al. 



19911 ). interpreted as indicative of the presence of a dense equatorial torus of dust. The opti- 
cal images show a central ionized region about 32" in diameter. CO maps show the presence 
of a collimated high velocity outflow in the polar direction and wit h [CO]/[H2]~ 5 x 10 -6 and 



a total molecular mass in the range 0.16-3.2 M Q ( ISahai et al.lll99ll ). Hubble Space Telescope 
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(HST) WFPC2 images in [N II], Ha and [O III] have revealed the existence of an intricate 



( 119971 ) found no evidence of radial variation for He, 0, N, Ne, and Ar. I Cox et al. Jl992h 



as 


Corradi et al. 


Cox et al. ( 


1992|) 



have detected several C-rich features at mm wavelengths, such as CN, HCO + , HCN and 
HNC, which ind i cate t he nebula is C-rich, although a C/O ratio of 0.6 is reported by 
Cohen fc Bartow! J2005h . 



IC 4406 is a relatively low electron density nebula. Values in the 400-2000 cm 3 range 



have been estimated using several different optical and i nfrared lines, wi t 
[S II] and [O III] doublets matching around 540 cm -3 (ILiu et al.l 12001 



1 values derived b y 



Wang et a]J|20QJ ). 



Its central star has a He II Za nstra temperature of 96800 K flPhillipd 120031 ) and its distance 
is probably ar ound 1.6 kpc flSahai et al.l Il99ll ) , although some authors claim it may be 
overestimated Jo'Dell et al.ll2002h . 



Gruenwald et al.l (119971 ) have modeled IC 4406 with a 3-D photoionization code and fit 
many observed line intensities assuming there is a torus around the central star. They find 
as a best fit a central star temperature of 8 x 10 4 K, luminosity of 400 L , torus density 
1500 cm -3 and nebular density 100 cm -3 . 

In general, comparisons of IR images of planetary nebulae, which trace the molecular gas 
and warm dust emission, to op tical line images, w hich trace the ionized gas, have shown the 
presence of similar structures (ILatter et al.lll995l ). leading to the conclusion that molecular 
and ionized gas spatially coexist in planetary nebulae, as well as dust grains, despite the 
different physical conditions these components are presumed to survive in. We have observed 
IC 4406 in the radio range to inspect the distribution of the ionized gas in its envelope and 
in the infrared to check for emission from the equatorial dust and molecular gas. 

In §2] we explain how we performed our observations and reduced the data; in $3] we 
show our results and in particular in §3.11 how we modeled the emission in the radio and 
infrared ranges; §3.21 compares our model results to the nebular parameter values obtained 
directly from the observational data; in §4] we summarize the present work. 
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2. 



Observations and Data reduction 



2.1. 



Radio observations 



Radio observations were 



performed at the Australia Telescope Compact Array 




on November 24, 2005 (17:00:00-08:00:00 UT) and December 11, 2005 (15:30:00-02:00:00 
UT), simultaneously at 4.8 and 8.6 GHz. The November run was performed with the array 
in 1.5C configuration, while for the December one the configuration was 6. OA. The adopted 
configurations are both linear but with different antenna positions giving maximum baselines 
of 4500 m (1.5C) and 5939 m (6. OA), minimum baselines of 77 m (1.5C) and 337 m (6. OA). 
The pre-calibration of the array was performed observing 0823-500, while the absolute flux 
calibrator was 1934-638. Another target was also observed during our two runs and the total 
on-target time was about 7 hours for each of the two. The phase calibrator chosen for IC 4406 
was 1431-48, which is 4°. 76 away from the target. The data were reduced with MIRIAD, 
following a standard reduction procedure as recommended in the MIRIAD User's Guide. 
The data from the two runs were combined into one dataset, obtaining a uv coverage from 
0.9 to 96 kA at 4.8 GHz and from 1.5 to 172 kA at 8.6 GHz. These correspond to an angular 
resolution of 2". 2 at 4.8 GHz and 1".2 at 8.6 GHz, while the largest observable structures 
(Largest Angular Scale) are 230" and 140" respectively. Such a setup suits the need to collect 
all the radiation from the target, whose maximum size, as previously mentioned, is about 



After combining the datasets from the two runs, we created a dirty map from the uv- 
dataset with the task INVERT. The image size was set to 1024 pixel in X band (8.6 GHz) 
but a larger size (2048 pixel) was necessary for the C band (4.8 GHz), to include a secondary 
source in the field and allow subsequent proper cleaning of the map. The cell size was set 
to 0".7 at 4.8 and 0".4 at 8.6 GHz, so that the beam would be extended over ~3 cells. 
We set the INVERT parameter options to mfs, double, so that the map could be obtained 
from a multi-frequency dataset without averaging in frequency and the size of the beam, 
created along with the map, would be the double of the map size, for a better performance 
of the cleaning algorithm in the following reduction step. To combine sensitivity to extended 
emission with sidelobe suppression, we used Brigg's weights, leaving the parameter sup unset 
and robust=l. Smaller values of the ROBUST parameter determine in our maps a worse 
signal-to-noise ratio. This setting gave a beam of 4".122x3".089 at 4.8 and 2".669x2".060 
at 8.6 GHz. 



1 The Australia Telescope Compact Array is part of the Australia Telescope which is funded by the 
Commonwealth of Australia for operation as a National Facility managed by CSIRO. 



100". 
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We then used the task CLEAN to deconvolve the dirty map from the synthetic beam. 
We set a gain of 0.1 and minpatch 257 for each CLEAN cycle, letting the task choose the 
proper algorithm, which at both frequencies was the CLARK one. We performed 3000 
CLEAN iterations, to reach the theoretical rms noise, as given in the output by INVERT. 

As at both frequencies our target appears to be quite resolved, the estimate of its flux 
has been performed on naturally weighted maps obtained using tapering (task INVERT: 
fwhm=15, for both frequency bands). This results in larger weights for the visibility points 
corresponding to short baselines, then limits the chance to miss extended flux in the final 
map, although possibly producing higher noise. Since the error in our measurement is 
primarily determined by the error in the absolute flux calibration, tapering is not an issue. 
With such procedure we measure at 4.8 GHz 103.3 ± 0.3 mJy and 92.5 ± 0.4 mJy at 8.6 
GHz, which, taking into account a typical 5% error in the absolute calibration, gives us the 
following final measurements: 103 ± 5 mJy at 4.8 GHz and 92 ±_5_mJy at 8.6 GH^. These 



values agree with the previous measurements in iMilne fc Allerl (119751 ) , performed with single 



dish telescopes, and therefore we can conclude that no flux is missing. 



2.2. Infrared observations 



Infrared observations were performed with the InfraRed Array Camera (IRAC) (IFazio et al. 



20041 ) onboard the Spitzer Space Telescope (Spitzer) at 3.6, 4.5, 5.8 and 8.0 /xm on March 06, 
2004 (UT 09:54:16.311) as part of the GTO program "Studying stellar ejecta on the large 
scale with SIRTF-IRAC" (AOR ID: 4414208). Six High Dynamic Range 30 sec dithered 
frames were obtained at each wavelength, for a total exposure time of 180 sec per channel. 

Basic Calibrated Data (BCD) were retrieved from Spitzer archive, clean ed to correct 



such artifacts as mux-bleeding and banding and then coadded using IRACproc ( jSchuster et al 



20061 ). For flux measurement the images were first converted from MJy/sr (IRAC BCD files 
are in units of MJy/sr) into Jy units using IRACproc. Then in each image four areas were 
boxed with a polygon using the task CGCURS in MIRIAD, the mean emission in each 
polygon was calculated then averaged to obtain an estimate of the background. The stan- 
dard deviation of the background mean value was assumed as an error. The whole nebula 
in each image was then boxed with a polygon of approximately 135" x 65". The size of 
the emitting nebula has been determined as the contour at the background+3<r level and 
is approximately: at 3.6 /xm 110" x 45", at 4.5 /xm 120" x 50", at 5.8 /xm 120" x 60" and 



2 The final error has been estimated as a = ^Jrms 2 + ((J ca iF) 2 , where a ca i is the 5% relative error in 
absolute calibration and F is the measured flux density. 
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at 8.0 jim. 125" x 65". The flux within each 135" x 65" polygon was summed up and the 
same procedure was adopted for the field star observed West of the central core, so that 
its flux was subtracted to the overall nebula's one. The result was corrected for extended 
emission according to the IRAC Data Handbook for aperture photometry. The error in the 
background estimate was assumed as the error of the flux measurement. We have measured 
at 3.6, 4.5, 5.8, and 8.0 /im the following fluxes respectively: 64 ± 2, 101 ± 4, 122 ± 7, and 
248 ± 8 mJy 

Table [1] summarizes the results of our radio and infrared observations. 



2.3. HST archive observations 

We have retrieved the calibrated Ha WFPC2 data from the HST web archive. The 
only image processing needed was to remove cosmic rays with the IRAF task CRREJ. The 
WFPC2 array is made of four different arrays, the central one (PC) having a pixel size of 
0".0455 and the other three arrays (WF) 0".0996. We selected the image obtained with 
the PC array and scaled it to match our radio beam (2".669x2".060) with the IRAF task 
GAUSS (Figj5]). The astrometry of the WFPC2 image was then modified with the IDL task 
HASTROM so that it could match our A TCA maps. HST images have a typical position 



uncertainty of ~0".5 (ILee fc Kwokl 120051 ) : we shifted the Ha image, to match the radio 



emission region, of 0".13 E and 0".2 S. 



3. Results and Analysis 

Our radio maps (Fig{T|) show the presence of a 42" x 56" (3cr level) emitting region at 
4.8 GHz, elongated in E-W direction; at a 10% of the peak level the size of the emitting 
region is restricted to 32" x 32". At 8.6 GHz the 10% of the peak level gives a size of 36" x 
40", in agreement with its 3<r level size. The maps do not show any N-S blobs of emission. 
What is seen is a clumpy emitting region that resembles what is observed by HST in Ha 



flO'Dell et al.ll2002h . 



Using the fluxes that we estimate at the two frequencies, we can calculate a spectral 
index a = —0.19 ±0.09, which matches the expected value of -0.1 for an optically thin radio 
shell. 

Fig|2] shows the IRAC images of the nebula. To properly view the central equatorial 
area, the images were plotted with a logarithmic scale having the peak flux and 3<r of the 
image as thresholds. Channel 1 and especially Channel 3 resemble the H2 emission image 
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in IStoreyl (119841 ) and this is probably because the v = 0-0 S(7) line at 5.5 1 nm falls withi n 
the 5.8 /iin band and several H 2 lines may contribute to the 3.6 /im image (IHora et a.ll2004l ). 
Channel 4 clearly shows the emission from the torus of dust surrounding the central star. 
Its size is about 28" x 20", elongated in N-S direction and the angular distance between its 
peaks is about 14" (the north peak is found at 14 /l 22 m 26 s .05, -44°08'54".64 and the south 
one at 14 ft 22 m 26 s .31, -44°09'08".14, with a position angle of 168°. 7). The overall size of 
the torus matches the approximate size of the nebula in the N-S direc tion (~ 3 0"), a lthough 



its peaks are much closer to the center than the H 2 blobs reported by IStoreyl (119841 ). whose 
separation can be roughly estimated as ~25". This indicates that the torus is partly shielding 
the molecular gas from the UV radiation from the central star. 

FigJH] is a combination of IRAC channels plotted with linear scale. Despite the lower 
resolution compared to Hubble images, IRAC is able to detect the faint emission from the 
neutral components in the envelope and reveals the structure of the elongated lobes. The 
IRAC images show filaments at different distances and inclinations from the central star, 
connected to the mass loss history of the nebula. The filaments that are closer in projected 
separation to the central star show relatively stronger blue (3.6 /im) emission, which may 
imply a higher temperature, being intrinsically closer to the central object. The overall 
structure observed in the envelope corresponds to the assumption that the central torus is 
the main collimating agent, as confirmed by the superposition of the 8.6 GHz and 8 fim 
images in FigJH 

T he scaled WFPC2 i mages match reasonably well our radio maps at 6 and 3 cm. Fol- 
lowing |Lee_^i£wokJ (120051 ) we also calculated the expected Ha image from our 8.6 GHz radio 
map and then an optical depth map (FigjSJ). The optical depth map confirms the clumpy 
nature of the central region and it also points out to a larger absorption toward the very 
core of the region, implying that dust can be present even in "close proximity" to the central 
star. 



3.1. The Spectral Energy Distribution of IC 4406 

To inspect the dust properties of our target we have collected literature data that, along 



with our observations, enable us to build the SED. We have retrieved 2MASS (ISkrutskie et al 



20061 ) images (J, H and Ks) from the 2MASS archive to measure our target's flux in such 
bands: the values in the Point Source Catalog neglect the emission from the extended en- 
velope. To estimate the flux and its error a procedure analogous to that used for IRAC 
data was applied. In each image four areas were boxed with a polygon, then the mean 
emission in each polygon was calculated, averaged to obtain the background and its stan- 
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dard deviation was taken as the flux density error for the selected image. The nebula itself 
was boxed with a polygon around its background+3cr contour, the flux within the polygon 
was calculated then background subtracted (the polygon size was the same as for the IRAC 
data). 2MASS fluxes were then con verted into UK I RT s ystem J, H and K magnitudes, 
so that color correction according to ISchlegel et al.l (119981 ) could be performed, adopting 
E(B-V)=0.19 (jGathier fc Pottaschlll988l ). The IRAS color corrected data were taken from 
Sahai et al.l (119911) . In our modeling we have assumed a central star temperature of 9 6800 
K ( jPhillipsI 120031 ). a distance of 1.6 kpc, and a luminosity of 170 L Q (ISahai et al.lll99ll ). 



We have modeled the SED separately for the radio (ATCA and single dish literature 
data) and infrared (IRAS, IRAC, 2MASS) emission. The radio data collected from the 
literature are listed in Table [2] with their references. For the former we have solved the 
radiation transfer equation in a spherical shell. The density distribution profile in the shell 
was determined fitting the density profile found in ICorradi et al.l (119971 ) , which gave us a 
radial variation as r _L252 . We introduced this density radial variation in our model shell, 
and found as a best fit to the data a density at the inner radius of 730 cm -3 , inner radius 
18", outer radius 24", having fixed a distance of 1.6 kpc and electron temperature 10 K. 
This gives us an optically thin spectrum down to 800 MHz, which confirms the status of 
our target as a fairly evolved object, and it also provides us with an estimate of the ionized 
gas mass of about 0.22 M . To estimate the free-free contribution from radio to near-IR 



wavele ngths, in our model we have calculated the Gaunt factor according to lKarzas fc Latter 



( 1196lh 



For the infrared range we have used the code DUSTY (llvezic et al.lll999l ) to solve the 
radiation transfer, assuming once again the nebula to be spherical. This simple assumption 
can provide reasonable constraints to the main properties of the envelope, since asymmetries 
in the density distribution would m ostly affect the optical part o f the SED, not the mid- 
IR/FIR region that we are modeling (jSanchez Contreras et al.ll20071 ). DUSTY does not allow 
the simultaneous treatment of more than one shell, yet its output can be used as an input 
in a second run, thus mimicing the shell structure of the nebula. 

One constraint to our model is the optical depth at a specified wavelength, which we 
can calculate as 



r„ = 2.03 x 10 1 



1 



e*B v (T d ) 

where T d is the dust temperature, 9 the angular radius of the nebula in arcsec, F v the 
flux density at the frequency v in erg cm " 2 s~ l Hz -1 , B^jT^) the Planck function at the 
temperature T d in erg cm -2 s _1 Hz -1 sr -1 (jGathier et al.lll986l ). 



The first attempts to fit the data were done with a standard MRN (jMathis et al.lll9771 ) 
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size distribution of the grains, with a m i n = 0.005, a max = 0.26, n(a) oc a q , q = —3.5 , density 
distribution in the shell as r -2 , r being the shell radius, a chemical composition of amorphous 
Carbon only and 0.1 as optical depth at 0.55 /im. The choice of an am-C only chemistry is 
due to the detection of several C-rich features mentioned in Section 1. 

By a first inspection of the observed data points, it was evident that the data could not be 
matched by a single dust component. In order to fit several components, we have performed 
our fit in steps, fitting first the lower wavelength data in a DUSTY run reproducing a hot 
inner shell, then giving the output of this run as an input to a second run of DUSTY. We have 
been able to reproduce the observed data assuming the dust envelope is made up of three 
shells, containing hot, warm and cold dust. The temperature of the cold component thus 
obtained was used in EqfTJto calculate r at 60 /jm. We chose this wavelength because at this 
wavelength all the flux seems to be due to one emitting component and cirrus contribution 
is negligible, which is not necessarily true at 100 //m. Having calculated the optical depth at 
60 with EqJT] we checked if this value matched the one given by DUSTY in its output. 
This was not the case. Then we used the tqo estimate as an input optical depth in DUSTY. 
This led to mismatch all the longer wavelength data points. We started changing the cold 
dust temperature, looking for a combination of T d and r 60 that would allow to match the 
data. We found that it was not possible to reach such a match with the specified set of 
input parameters, the DUSTY fluxes at larger wavelengths being larger than the observed 
ones. We have then tried to change the density distribution dependence on the radius: such 
exponents as -3, -1, -0.5 were tested but none resulted into a good match to the observations. 
Our second attempt was changing the grain size: we noticed that the presence of larger grains 
in the model (up to 6.5 jim) could effectively modify the reproduced data. 

Since DUSTY makes use of spherical geometry, we have assumed an effective radius of 
45", corresponding to the radius of a circle having the same area as the dust ellipse observed 
in our IRAC images. The shell relative thickness parameter in DUSTY has been calculated 
to reproduce this angular size. The final set of parameters for our best fit is reported in 
Table El 

FigJT] shows the fit components to the observational data points and the combination of 
the infrared and radio fits. 

We notice that, if in our free-free model we had used the usual radio approximation 
of the Gaunt factor, the model would predict around 2 /im a lower level of emission than 
observed (as can be seen in Figure El where the single emitting components are plotted), 
which might be interpreted as due to a fourth hotter component of dust missing in the 
model. Our proper estimate of the Gaunt factor shows how the free-free contribution in the 
near-IR is actually non negligible and allows us to achieve a good fit to the data points in 
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this range. The fairly larger flux measured in K band can be explained when considering 
that H 2 emission has been detected in IC 4406 (jStoreylll984j ) and se veral lines may fall within 
the K band filters, along with ionized g as lines such as Br7. In fact iPhillips fc Ramos-Larios 
( 120051 ) and iRamos-Larios et al.l (120061 ) have explained the excess in Ks band in terms of H 2 
e mission, i n par ticular the latter show how the Ks band 2MASS image matches the H2 image 



111 



Storev! (I1984h . 



3.2. Nebular parameters 

Our models enable us to estimate such nebular parameters as electron density, ionized 
gas mass and dust ma ss. The dust mass can be calculated from the DUSTY output following 



Sarkar fc Sahail (120061 ) as 



M d = 4nR 2 Y^- (2) 

&100 

where R is the inner radius of the emitting shell in cm, Y is the thickness of the shell relative 
to R, T100 and ki o the optical depth and absorption coefficient at 100 //m. We can now use 
our DUSTY output for R and rioo, which, for the 57 K more external shell, are estimated a s 



3.6xl0 16 cm and 2.34xl0~ 5 ; Y is 30 and k 100 =92 cm 2 g" 1 , calculated following |jura| (119861 ). 
We thus obtain a dust mass of 6xl0 -5 M . Our radio model gives instead an average 
electron density of 450 cm -3 and ionized mass of 0.29 M , as mentioned in the previous 
section. 

These values can be compared to those derived by equations that directly use the ob- 
served fluxes. From our radio observations it is possible to derive the Hp flux and electron 
density, which can then be used to estimate the ionized mass of the nebula. 

jj S^.sGHz -2 -1 ,n\ 

Hr — _ „ cor ; r 7 era cm s (3) 

p 2.83 x 10 9 t°- 53 [l + (1 -x)y + 3.7xy] y KJ 



+0.88 u 

n e = 2.74x10^—^ cm- 3 (4) 

11.06 d 2 t°- 88 TT „ 

M ion = Hp M 5 

n e 

In EqjniH^ is the FL3 line flux in units of 10 -11 erg cm -2 s _1 , n e the electron density in 
cm -3 , t the electron gas temperature in units of 10 4 K, d the distance to the star in kpc. In 
EqJH Hg is like in the previous equation, 9 is the ionized gas radius as deduced from the 4.8 
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GHz radio map (we have used a value of 24") in arcsec, e is the filling factor, for which we 
have used an average value of 0.6, t as before. In Eqj3], S^sghz is the 4.8 GHz flux density in 
J y, t as before, x is He ++ /He = 0.121, y is He /H = 0.1 3 2, cal culated from the abundances 



in 



Corradi et all (119971 ). EqEl and Hare from lPottaschl (11984 1 We find 0.21 M , 418 cm 



and 3.09x10 11 erg cm 2 s 1 for the ionized mass, electron density and Hg flux respectively, 
which closely match our radio model's results. 



The dust mass can be estimated according to iPottasch et aJj (119841 ) as 

4 ap d 2 F y 



M. 



dust „ 



3Q v B(u,T t 



9 



(6) 



dust J 



Following lJural (119861 ). we consider 1.1 as a representative value of the power- law dis- 
tribution of the emissivity of carbon grains in the infrared. We can thus calculate the dust 
emissivity at 60 /im assuming average grain radius a = 10~ 5 cm and density p = 3 g cm -3 , 
which results into Qqo = 1-107 x 10~ 3 . Considering the flux at 60 pm as due only to the cold 
component at 57 K in our best fit, we calculate a dust mass of about 2.8 x 10 -4 M Q . We 
also calculated the dust mass of the other emitting components in our fit, but they resulted 
to be negligible when compared to the cold dust: we find for the 700 K component 5 x 10~ n 
M Q (calculated using the 4.5 pm flux) and for the 200 K component 4 x 10 -8 M (using 
the 12 pm flux). The dust to gas mass ratio can be estimated as Md/M g = 1.3 x 10~ 3 . See 
Table H] for a summary of the parameters we have derived. 

We notice that the cold dust mass value we estimate by DUSTY is one order of mag- 
nitude smaller than that calculated with Eqj6j This could perhaps be due to the approx- 
imations intrinsic to both the DUSTY modeling and the derivation of Eqj6] (i.e., spheri- 
cal shape, singl e comp onent chemistry, physical knowledge of dust opacity). For example, 
Sarkar fc Sahail ( 20061 ) noticed how in general the DUSTY SEDs are not very sensitive to 
cooler dust at large radii, as demonstrated by the fact that large differences in shell relative 
thickness values do not determine drastically different SEDs. Therefore it is possible that 
the code itself is underestimating the amount of dust. 

Another issue concerning the use of DUSTY is the assumed spherical symmetry. As pre- 
viously mentioned, one-dimensional models have been utilized in the literature for modeling 
the mid- and far-IR SEDs in planetary nebulae, since asymmetrie s in the distribution and 
orien tation have less of an effect than in the optical part of the SED ( jSanchez Contreras et al. 
2007h . 



Using more realistic geometries with other codes would be expected to provide more 
accurate results, but assumptions must be usually made on the 3D geometry. The advan- 
tages of a better match to the geometry of the nebula given by these codes can be heavily 
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affected by such assumptions, so a careful exploration of parameter space consistent with 
the observations must be performed in order to guide and interpret the modeling results. 

The use of one-dimensional codes as DUSTY is therefore still valid in determining the 
mean properties of an envelope, within the intrinsic errors of any modeling process. 

4. Summary 

We have observed IC 4406 in the cm and 3-10 micron ranges. Our radio observations 
have confirmed the presence of the complicated maze of lanes already observed in Ha in the 
central region of the nebula and have enabled us to calculate several nebular parameters, 
whose values match the classification for this target as an evolved planetary nebula, in 
particular its low dust to gas mass ratio and density. IRAC imaging has revealed the presence 
of filaments in the nebula that were not detected in previous observations. 

Our IRAC measurements, combined with literature data at longer and shorter wave- 
lengths, have enabled us to study the SED of the PN IC 4406 and reproduce it with DUSTY. 
This has revealed that three different dust components are needed to model the data, with 
temperatures ranging from 57 to 700 K. It has also been necessary to include in the model 
slightly larger grains than in the standard MRN composition (up to 6.5 /im) to account for 
the calculated 60 /im optical depth. The main limits of our modeled curve are the spherical 
geometry assumed in DUSTY and the lack of data in the mm and sub-mm ranges, which 
would give a constraint on the slope of the curve. As we have observed during our trials 
with DUSTY, the slope of the SED in the sub-mm range changes with the maximum size 
of the grains included in the model. Unfortunately, in this range observations are available 
only for a few stars so far: none for our target. 

We can speculate that in such a diversified dust environment, as we find in IC 4406, 
further lower temperature components may exist and future high sensitivity, high angular 
resolution observations will give a fundamental contribution to understand the physics of 
circumstellar envelopes in planetary nebulae. 
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Facilities: ATCA, Spitzer (IRAC) 



Table 1: IC 4406 measured fluxes. 



Spitzer 


3.6 fim 


4-5 \xm 


5.8 \xm 


8.0 fim 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


IC 4406 


64 ± 2 


101 ± 4 


122 ± 7 


248 ± 8 


ATCA 




4.8 GHz 


8.6 GHz 








(mJy) 


(mJy) 




IC 4406 




103.3 ± 0.3 


92.5 ± 0.4 




1431-48 




1020 ± 1 


650 ± 1 





13 



Table 2: IC 4406 radio fluxes from the literature. 



v (GHz) 


S„ (mJy) 


0.843 


104 ± 3 a 


2.7 


150 ± 60 b 


5 


110 ± 15 c 


14.7 


84 ± 8 d 



"Mauch et al. ('< 


>003 


) 


1 Milne & Webster ( 


1979 N 


c Milne & Aller 


(1975) 


c Milne & Aller 


(1982) 



Table 3: DUSTY fit parameters. 

Parameters Values 

Chemistry 100% amorphous Carbon 

Central source Blackbody at 96800 K 

Density distribution oc r -2 

Grain size distribution oc r -35 , a min = 0.005 juiii, a max = 6.5 /im 

Hot component T d = 700 K, r 60 = 4.5 x 10" 6 , R in = 7.17 x 10 13 cm, Y=14790 

Warm component T d = 200 K, r 60 = 1.0 x 10" 5 , R in = 1.40 x 10 15 cm, Y=750 

Cold component T d = 57 K, r 60 = 6.8 x 10" 5 , R in = 3.62 x 10 16 cm, Y=30 



Note. — Y is the shell thickness relative to the inner radius, which is determined by the code as the layer 
having temperature Td] Ri n is the distance to the central star of the layer having temperature T d ; tqq is the 
optical depth at 60 /im. 



Table 4: Calculated nebular parameters. 



Empirical Models 
M ion (M ) 0.21 0.22 
M dust (10- 4 M ) 2.8 0.6 
M dust /M wn (10- 4 ) 13 3 
Gas density (cm" 3 ) 418 730 



Note. — As "Empirical" we mean values calculated with the equations in and with "Models" those 
deduced from the modeling explained in §3.11 The model gas density is the value at the inner radius. 
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Fig. 1. — Radio maps of IC 4406 (top: 4.8 GHz, bottom: 8.6 GHz). The contours indicate 
the (-3, 3, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110)<r and the (-3, 3, 6, 9, 12, 15, 18, 21, 
24, 27) a levels for the top and bottom map respectively, with a ~0.05 mJy for both maps. 
Both maps were produced with uniform weighting and ROBUST=l. The synthetic beam is 
shown in the bottom left of each map and the flux density unit is Jy/beam. 
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Fig. 2. — The emission observed with IRAC. The scale is logarithmic: the maximum (white) 
is the emission peak in the area, the minimum (purple) is 3<r (maxima are respectively for 
3.6, 4.5, 5.8 and 8.0 /im: 3.97, 5.85, 6.75 and 15.07 MJy/sr, while a is 0.09, 0.14, 0.14, and 
0.52 MJy/sr). Clockwise from top left: 3.6, 4.5, 5.8, 8.0 /im. 
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Fig. 3. — Images obtained with IRAC. Top: overlay of IRAC channels (blue is 3.6, green 
4.5, orange 5.8 and red is 8.0 /mi). Middle: as in the top image but not including 5.8 //m 
channel, for better viewing the central core emission. Bottom: 5.8 /mi image scaled for a 
better view of the filament structures, some of which are pointed by the arrows. 
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RIGHT ASCENSION (J2000) 

Fig. 4. — IRAC channel 4 image is shown in false color (units are MJy/sr) and overlayed are 
the 8.6 GHz contours as in Fig{TJ 



-18- 




Fig. 5. — HST WFPC2 image scaled to match our 4.8 (top) and 8,6 GHz (bottom) maps. 
The units are DN and the absolute values of the levels -3, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 
33, 36, 39, 42. 
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Fig. 6. — Optical d epth map of the cen tral region, obtained from the 8.6 GHz map and the 
Ha image following iLee fc Kwokl (120051 ). 
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Fig. 7. — SED of IC 4406 from near-IR to radio wavelengths. The best fit is obtained 
summing the DUSTY and free-free continuum model outputs, and is shown as a thick solid 
line. The single fit components are also shown. The dot curve is the result of a model with 
the central star and the 700 K component, the short dash curve includes the central star, 
700 K and 200 K components and the dot- dash curve includes the central star and the 700, 
200, and 57 K dust components, the long dash curve is the free-free continuum. Error bars 
are shown for all data, some are smaller than the data point. 
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